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Abstract 



Global and local regularities of functions are analyzed in anisotropic function spaces, 
under a common framework, that of hyperbolic wavelet bases. Local and directional 
I regularity features are characterized by means of global quantities constructed upon the 

■ coefficients of hyperbolic wavelet decompositions. A multifractal analysis is introduced, 

that jointly accounts for scale invariance and anisotropy. Its properties are studied in 
depth. 
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>: 

^ ■ 1 Introduction 

■ Natural images often display various forms of anisotropy. For a wide range of applications, 

CN ■ anisotropy has been quantified through regularity characteristics and features that strongly 

differ when measured in different directions. This is, for instance, the case in medical imaging 
(osteoporosis, muscular tissues, mammographies,...), cf. e.g. [15i|l4J, hydrology ^44J) fracture 
^ ■ surfaces analysis [22] v • • • For such images, a key issue consists first in describing, within a 

^ ! suitable framework, the anisotropy of the texture, and then in defining regularity anisotropy 

parameters that can actually and efficiently be measured via numerical procedures and fur- 
ther involved into e.g., classification schemes. This requires the design of a mathematical 
framework that allows to define and estimate these parameters. Such a program can be split 
into several questions, some of them having already been either solved or, at least, patially 
addressed. 

A first issue lies in introducing global and local notions of anisotropic regularity, which 
emcompass and extend (isotropic) regularity spaces, such as Sobolev or Besov spaces, and 
the classical notion of pointwise Holder regularity. To model anisotropy, the particular set- 
ting of an anisotropic self-similar field driven by two parameters (an anisotropy matrix and a 
self-similarity index) has been introduced and studied in [14J, where it is used as a relevant 
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model to describe osteoporosis. In [21J, the question of defining in a proper way the concept 
of anisotropy of an image in relation to its global regularity has been addressed. It has no- 
tably been shown that these two parameters can be recovered without a-priori knowledge of 
the characteristics of the model, by studying the global smoothness properties of the process. 
Furthermore, some of the properties characterizing anisotropy are revealed by the regular- 
ity of the sample paths when analyzed with functional spaces well-adapted to anisotropy: 
Anisotropic Besov spaces. 

This preliminary study thus showed the central role that such spaces should play in the 
mathematical modeling of random anisotropic textures. The introduction of these spaces 
traces back to the study of some PDEs, cf. e.g. [51], for the study of semi-elliptic pseudo- 
differential operators whose symbols have different degrees of smoothness along different di- 
rections, the reader is also referred to [2], and references therein, for a recent use of such spaces 
for optimal regularity results for the heat equation. Other types of directional function spaces 
have also been considered, cf. e.g. [I6] for the variant supplied by Hardy directional spaces. 

A second crucial issue consists in obtaining a simple characterization of these spaces on 
a "dictionary". Indeed, the challenge is to measure the critical exponent of any image in 
anisotropic Besov spaces for different anisotropics using the same analysis tool. Wavelet 
analysis is well-known to be an efficient tool for measuring smoothness in a large range of 
functional spaces (cf. [39] for details). Here, however, the main point is that the anisotropy of 
the analyzing spaces must not be set a priori to a fixed value but instead be allowed to vary. 
Specific bases are thus looked for, which would serve as a common dictionary for anisotropic 
Besov spaces with different anisotropics. It is natural that one should use some form of 
anisotropic wavelets such as curvelets, bandelets, contourlets, shearlets, ridgelets, or wedgelets 
(see e.g. [34J for an thorough review of these representation systems and a comparison of their 
properties for image processing); natural criteria of choice being, on the mathematical side, 
that these variations on isotropic wavelets supply bases for the corresponding anisotropic 
spaces, and, on the applied side, that practically tractable procedure can be devised and 
implemented to permit the characterization of real-world data according to these function 
spaces. 

Many authors addressed this problem, and proposed different solutions, depending on the 
precise definition of anisotropic space they started with, as well as on the anisotropic basis 
they used, cf. e.g. [231 1331 ISSl [Ml [52] and also the recent papers [28j by G. Garrigos and 
A. Tabacco, and |31j by D. Haroske and E. Tamasi, which contain numerous references on 
the subject. Note also that, in several particular type of anisotropy was considered: 

Parabolic anisotropy (where a contraction by A in one direction is associated with a contrac- 
tion by along the orthogonal direction, [371 [l6l H2] . and references therein, in particular 
for applications to directional regularity), the corresponding dictionaries being in that case 
curvelets or contourlets (corresponding to the Hart-Smith decomposition in the continuous 
setting, see |49]). This particular choice of anisotropy was motivated from application to 
PDEs (see e.g., [29l [20] where curvelets and ridgelets are used for the study of Fourier inte- 
gral operators, with applications to the wave equation) but is no longer justified when dealing 
with images, where no particular form of anisotropy can be postulated a priori. To the op- 
posite, figuring out the precise form of anisotropy present in data is part of the issue. This 
argument also implies that one should not restrict analysis to tools that match one specific 
type of anisotropy, but rather that to tools embracing all of them simultaneously, in order to 
be able to detect that that suits data. 

From now on, two possible solutions for this problem will be focused on: 
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• One is supplied by anisotropic Triebel bases, see [52], that are constructed from the stan- 
dard wavelet case through a multiresolution procedure, tailored to a specific anisotropy. 
The collection of these bases does not constitute a frame. However, for a fixed anisotropy, 
simple characterizations of anisotropic Besov spaces have been supplied within this sys- 
tem. Such characterizations can thus be used as a building step to construct a multi- 
fractal formalism [?]. This is further detailed in Section [3l 

Triebel bases provide a powerful tool to deduce results on anisotropic Besov spaces, for 
a fixed anisotropy. In particular, it enables to show that these spaces are isomorphic to 
the corresponding isotropic Besov spaces. Further, some results such as embeddings or 
profiles of Besov characteristics can be obtained, via the transference method proposed 
by H. Triebel. However, when it comes to understand the link between different forms 
of anisotropy - in term of function spaces by example - this tool remains of limited 
interest. Indeed, the knowledge of the expansion of a function in one basis gives a priori 
neither information about its expansion in an other basis nor about its belonging to all 
anisotropic Besov spaces. 

• Another possible decomposition system is supplied by hyperbolic wavelets, introduced 
in various settings under different denominations (standard, rectangular or hyperbolic 
wavelet analysis) notably in image coding (see [51]), numerical analysis (see jT3], [12]) 
and in [23j,[33j for the purpose of approximation theory. They are simply defined as 
tensor products of ID wavelets, yet allowing different dilations factors along different 
directions, as opposed to the classical discrete wavelet transform that relies on a single 
isotropic fixation factor. This key difference enables the study of anisotropy. Hyperbolic 
wavelet basis form a non-redundant system by construction, and contain all possible 
anisotropics. Hyperbolic wavelet bases have thus been used in statistics for the purpose 
of adaptive estimation of multidimensional curves. Notably, it has been proven in two 
seminal articles [U] and |10] that nonlinear thresholding of noisy hyperbolic wavelet 
coefficients leads to (near)-optimal minimax rates of convergence over a wide range of 
anisotropic smoothness classes. The reader is also referred to the recent work of F. 
Autin, G. Claeskens, J.M. Freyermuth [3] where this problem is considered from the 
maxiset point of view. Other interesting applications of hyperbolic analysis can also be 
founded in [5j , [4j , [8j , [6j and [7j where hyperbolic wavelet decompositions of Fractional 
Brownian Sheets and Linear Fractional Stable Sheets are given and are used to prove 
many sample paths properties of these random fields (smoothness properties, Hausdorff 
dimension of the graph). 

The key feature of hyperbolic wavelet bases is that they provide a common dictionary 
for anisotropic Besov spaces. This result is stated in Theorem 12.21 of Section [2j The 
critical exponent in anisotropic Besov spaces will be related to some norms of the 
hyperbolic wavelet coefficients. These mathematical results yield an efficient method 
for the detection of anisotropy, as detailed in a companion article, where numerical 
investigations are conducted, [45j . 

In the present article, it has been chosen to explore the possibilities supplied by the hyper- 
bolic wavelet transform to investigate directional regularity, both in global (anisotropic Besov 
spaces) and local (directional pointwise regularity) forms. The underlying motivation is to 
develop a multifractal formalism relating these two notions (just as the standard multifractal 
formalism relates the usual Besov spaces with the notion of (anisotropic) Holder pointwise 
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smoothness, see [35] and references therein). It also ahns at obtaining a numerically stable 
procedure that thus permits to extract the anisotropic features existing in natural images as 
well as information related to the size (fractional dimensions) of the corresponding geometrical 
sets. 

Before proceeding further, let us motivate the choice of hyperbolic wavelets against Triebel 
bases. For a fixed anisotropy, one can argue that Triebel bases display slightly better mathe- 
matical advantages: An exact characterization of anisotropic Besov spaces, as shown in [52] , 
and a characterization of pointwise smoothness as sharp as in the isotropic shown 
by H. Ben Braiek and M. Ben Slimane in [lOj. However, a first purpose of the present con- 
tribution is to show that these two important properties hold almost as well for hyperbolic 
wavelets: In Section [21 "almost characterizations" (i.e., necessary and sufficient conditions 
that differ by a logarithmic correction) of anisotropic Besov spaces are obtained. Further- 
more, if one is not only interested in analysis, but also in simulation, this slight disadvantage 
(a logarithmic loss, which in applications can not be detected) is overcompensated by the 
advantage of using a basis instead of an overcomplete system. Indeed, generating a random 
field with prescribed regularity properties requires the use of a basis (using an overcomplete 
system cannot guarantee a priori that the simulated field with coefficients of specific sizes has 
the expected properties, since nontrivial linear combinations of the building blocks may van- 
ish). A contrario, with the hyperbolic wavelet basis, one can easily provide toy examples with 
different multifractal spectra depending on the anisotropy. Our being jointly motivated by 
analysis and synthesis motivates the choice of a system that permits an interesting trade-off 
among directional wavelets, in terms of mathematical efficiency and numerical simplicity and 
robustness, both on the analysis and synthesis sides. The practical relevance of the mathe- 
matical tools introduced and studied here are assessed in a companion paper [35] . 

Let us now further compare Triebel and hyperbolic wavelet bases in terms of pointwise 
directional smoothness. First, note that this notion has been the subject of few investigations 
so far: To our knowledge, the natural definition which allows for a wavelet characterization 
was first introduced by M. Ben Slimane in the 90s, see [TTj , in order to investigate the mul- 
tifractal properties of anisotropic selfsimilar functions. Partial results when using parabolic 
basis (i.e., curvelets and Hart-Smith transform) have been obtained by J. Sampo and S. 
Sumetkijakan see [37^ H6| H2] and references therein. A generalization and implications in 
terms of sizes of coefficients on directional wavelets (the so-called "anisets" , which are a mix- 
ture of of the wavelet and Gabor transform, where the wavelets can be arbitrarily shrunk 
in certain directions) were also worked out in [36]. Finally, an "almost " characterization of 
pointwise directional regularity was recently obtained by H. Ben Braiek and M. Ben Slimane 
in [?] on the Triebel basis coefficients, where the basis is picked so that its anisotropy pa- 
rameter is fitted to the type of directional regularity considered. In Section [21 we will obtain 
a similar result, but relying on the coefficients of the hyperbolic wavelet basis, thus paving 
the way to the construction of a multifractal formalism. An important difference with [?] is 
that, here, a single basis fits all anisotropies. Therefore, as in the case of Besov spaces, the 
advantage is that no a priori needs to be assumed on the particular considered anisotropy. 
This thus can be used as a way to detect the specific anisotropy which exists in data at hand, 
rather than assuming a priori its particular form beforehand. Note that other decomposition 
systems have also been used for the detection of local singularities, see for instance [241 130] 
where shear lets and wedgelets are used for the detection of discontinuities along smooth edges. 
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Let us now come back to the anisotropic self-similar fields considered in |451 [2T] . Such 
exactly selfsimilar models are somewhat toy examples, and, though testing regularity indices 
on their realizations is an important validation step, their study could prove misleadingly 
simple (just as, in ID, fractional Brownian motion is too simple a model to fit the richness of 
situations met in real- world data). Natural images are indeed likely to consist of patchworks 
of different kinds of deformed pieces and therefore, can be expected to exhibit more complex 
scale invariance properties, and only in an approximate way. A natural setting to describe 
such properties, where different kinds of singularities are mixed up, is supplied by multifractal 
analysis. The next step is therefore to combine both anisotropy and multifractal analyses. To 
this end, a new form of multifractal analysis is introduced, based on the hyperbolic wavelet 
coefficients, and relating the global and local characterizations of regularity. It allows to take 
into account both scale invariance properties and local anisotropic features of an image. Thus, 
it provides a new tool for image classification, seen as a refinement of texture classification 
based on the usual isotropic multifractal analysis, as proposed for instance in [1],[35]. Section[3] 
is devoted to the introduction of this new framework: A new multifractal formalism, referred 
to as the hyperbolic multifractal formalism. It allows to relate local anisotropic regularity of 
the analyzed image to global quantities called hyperbolic structure functions as commonly 
done in multifractal analysis. Note that alternative multifractal analysis and multifractal 
formalism were introduced by H. Ben Braiek and M. Ben Slimane in [?], based on Triebel 
basis coefficients. In their approach, a particular anisotropy is picked, and the corresponding 
basis is used. As above, the main difference between our point of view and theirs is that 
we do not pick beforehand a particular anisotropy: Therefore, the approach proposed here 
does not rely on any a priori assumptions on data, and can thus be used when anisotropics 
of several types are simultaneously present in data. 

Finally, detailed proofs of all the results stated in Sections[2]and[3]are provided in SectionSl 

2 Anisotropic global regularity and hyperbolic wavelets 

We first focus on the measure of anisotropic global regularity using a common analyzing dic- 
tionary: hyperbolic wavelet bases. Here, we start by providing the reader with a brief account 
of the corresponding functional spaces. Thereafter, we recall some well-known facts about 
hyperbolic wavelet analysis (cf. Section \2.2\i . The main result of the present section con- 
sists of Theorem 12.2^ proven in Section 14. 1^ which allows to determine the critical directional 
Besov indices of data by regressions on log-log plot of quantities based on hyperbolic wavelet 
coefficients (see Section [2.21 for a precise statement). 

2.1 Anisotropic Besov spaces 

Anisotropic Besov spaces generalize classical (isotropic) Besov spaces, and many results con- 
cerning isotropic spaces have been extended in this setting, see [iHl [17] for a complete account 
on the results used in this section, and [Ml [53] for detailed overviews on anisotropic spaces. 
Note in particular that these spaces are invariant by smooth diffeomorphisms on each coor- 
dinate, an important requirement for image processing. 

Anisotropic Besov spaces verify (asymptotically in the limit of small scales) norm invari- 
ances with respect to anisotropic scaling, we, therefore, start by recalling this notion. Let 
a = (01,02) denote a fixed couple of parameters, with ai, 02 > and ai + 02 = 2. In the 
remainder, such couples will be referred to as admisible anisotropics. Such couples quantify 
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the degree of anisotropy of the space (ai = 02 = 1 corresponding to the isotropic case). For 
any t > and ^ = (^1,^2) £ I^^i we define anisotropic scahng by = (t^^^^i, ^"^^2)- Note 
that, in this definition and in the following, the coordinate axes are chosen as anisotropy di- 
rections. This particular choice can of course be modified by the introduction of an additional 
rotation (as envisaged e.g., in [35]). 

Anisotropic Besov spaces may be introduced using an anisotropic Littlewood Paley anal- 
ysis, which we now recall. Let (Pq>0 belong to the Schwartz class 5(M^) and be such 
that 



and 

For j € N, we define 
Then, 



ip^{x) = 1 if sup < 1 , 

i=l,2 



^^ix) = if sup |2-"»ei > 1 • 

't=l,2 



j=0 

and {^'j)j>o is called an anisotropic resolution of the unity. It satisfies 

supp((^^)ci??, suppiift) C R'^+ARj , 

where 

={e = (ei,e2)eM2; sup ie^i<2-'^}. 

i=l,2 

For / G cS'(M2) let 

A'^f = T-^ {^^y 

The sequence (A"/)j>o) is called an anisotropic Littlewood-Paley analysis of /. The anisotropic 
Besov spaces are then defined as follows (see [I8l [T7]). 

Definition 2.1 The Besov space B^''^^^^^^i^{R'^), forO<p< +(x>, < q < +00, s,/3 £R, is 
defined by 

log 1^(1^') = {/ e s'iR'y, r ^"2^1^'^ w}j < +00} . 

This definition does not depend on the resolution of the chosen unity ipQ and the quantity 

\j>0 

is a norm (resp., quasi-norm) on Bp^(R^) for 1 < p, q < +00 (resp., < p, q < 1). 
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As in the isotropic case, anisotropic Besov spaces encompass a large class of classical 
anisotropic functional spaces (see [53] for details). For example, whenp = q = 2 and (ai, 02) G 

is an admissible anisotropy, let us consider s > such that s/ai and s/a2 are both integers, 
then the anisotropic Sobolev space 

//^'"(]R2) = {/ e l2(m2) such that G L2(m2) and £ L^m^)} , 

coincides with the Besov space i?2'2 (^^)- 

In the special case where p = q = co, the spaces B^ooO^'^) are called anisotropic Holder 
spaces and are denoted C^{^^ |„ (M^). These spaces also admit a finite difference characterization 
that we now recall (see also |53j for details). 

Let (61,62) denote the canonical basis of M?. For a function / : — t- M, £ G {1)2} and 
t G M one defines 

Ayix) = fix + tee) - f{x) . 

The difference of order M, M > 2, of function /, along direction e^, is then iteratively defined 
as 

A*'{/(x) = A,,Aj'^-V(x) . 

One then has: 



Proposition 2.1 Let a = (ai,a2) G (M^) such that ai + 02 = 2, s > 0, u G M and 

-is,a , 
1log|"* 



/ : — 7- M. The function f belongs to the anisotropic Holder space Cf!°' |u (M^) if 



\We fix)\\L-°{K^) 
+ / sup < +00 , 

^t>o |ih/'^^|log(|t|)|« 
where for any £ G {1, 2}, Mi = [s/ai] + 1- 

2.2 Hyperbolic wavelet characterization of anisotropic Besov spaces 

We state our first main result which consists of an hyperbolic wavelet caracterization of 
anisotropic Besov spaces. 

We first need to recall the definition of the hyperbolic wavelet bases as tensorial products 
of two unidimensional wavelet bases (see [23] ) and second state Theorem 12.21 further proven 
in Section HTTl 

Definition 2.2 Let tp denote the unidimensional Meyer wavelet and ip the associated scaling 
function. The hyperbolic wavelet basis is defined as the system {V'ji,j2,fci,fc2' (^15^2) G (Z+ U 
{-1})2, {ki,k2) G where 

• if ji,j2 > 0, 

V'ji,j2,fci,fc2 (2:1,^:^2) = ip{2^''xi - ki)ip{2^^X2 - k2) . 

• if ji = -1 and j2 > 

ip-ij^MMi^^^^'^) = ^i^i ~ ki)ipi'^^'^X2 - k2) . 
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• ifji > andj2 = -1 

V'ii -i,fei,fe2(a^i)^2) = i'i2^''xi - ki)ip{x2 - k2) ■ 

• if h =32 = -1 

■0-1 -l,ifci,fe2(a^l>a;2) = '^{Xi - ki)(p{x2 - k2) . 

For any f G (S'(M^), one then defines its hyperbolic wavelet coefficients as follows: 

^h,32MM = 2^'^^^ < f,i^ji,j2M,k2 > ifhih > , 

Cji -i,fci,fc2 = 2-'' < f,ipjij2,kuk2 > if h > and 32 = -1 , 

C-i^j2MM = 2^^ < f,4'ji,j2MM > if h = -1 and j2 > , 

C-i-iMM = < /> V'iij2,fci,fc2 > ifh = h = -1 • 

Remark 2.1 We chose a -normalization for the wavelet coefficients, known to be well- 
matching scale invariance. 

The main result of this section is an hyperbohc wavelet characterization of the spaces 
j^gl^(R^), up to a logarithmic correction. In the sequel, some notations are needed. For 

any j = {ji^h) e (NU {-1})2, let us define: 

\/p 

\P 




Let a = {0.1,0.2) be an admissible anisotropy, one also defines the following subsets of 

rf ^)(a) = {(ji,j2) G n\ [{j - - 1 < ii < [jai] + 1 and < 32 < [{j - 1)02] - 1}, 

rJ.^^V) = iUuh) G < ii < [{j - l)ai] - 1 and [(j - 1)02] - 1 < J2 < [jaa] + 1}. 



j 

and 



r;.^^)(a) = {(ii,i2) e [(j - l)ai] - 1 < ii < [jai] + 1 and [(j - 1)02] - 1 < j2 < [ja2] + 1} , 



(1) r,- (a) = rf ^) (a) u rf ^) (a) u rf ^) («) . 

Let us now state our hyperbolic wavelet characterization of anisotropic Besov spaces: 

Theorem 2.2 Let a = (0:1,02) be an admissible anisotropy, {s,/3) G and{p,q) G (0, +00]^. 
LetfeS\R^). 

1. 5ei^(p,g) = max(l/p-l,0) + max(l-l/g,0). // 

. (Jl+J2)9 



jql3ip,q)-l3g2jsg ^ 2" 
f/ien / G 5*'^|j^g|^(R^) (^mi/i usual modifications when q = 00). 




2. Conversely, 



(a) Ifq<ooandfe B'^2\iog\^(^ ) ^^^n 

||c,,,,,||,^J <+oo. 



2 

VJGNo (ii,j2)er,(a) 



^^^^//e^;;::.,iiog|.(^') ^^^^ 

max 7 '^^2-'^ max 2 p lie,-, . .\\pp < +oo . 
jeNo-" 0-i,i2)er,(a) n Jij2, , 

Proof. Theorem 12.21 is proven in Section [4.H □ 

In particular, for p = q = oo, the following "almost" characterization of anisotropic Holder 
spaces by means of hyperbolic wavelets holds: 

Proposition 2.3 Let a = (01,02) o,n admissible anisotropy, S M? and f £ S'{M?). 

(i) If f £ C'^'"(M^) then there exists some C > such that for all j G N U { — 1} and any 

{31,32) e rj(a), 

(2) ||c,,,,,,,||,^ <C2-^'^ 

(a) Conversely, assume that there exists some C > such that for all j E NU{— 1} and any 
{31,32) e rj(a) 

(3) ll'^ii j2,-,-lk°° — , 
then f G C^'"(M2). 

In the special case where p = q = 2, that is if we consider anisotropic Sobolev spaces, there 
is no logarithmic correction: 

Theorem 2.4 Let a = (01,02) «^ admissible anisotropy, s G R. Let f £ S'(M?). The two 
following assertions are equivalent: 

(i) feH^'-{M.') = B-^l^{R^). 
(n) 

\ 1/2 

^2^- E 2-(-+-) ||c,,,,,||,0 <+oo. 
Wo (jij2)erj(a) / 

Proof. Theorem 12.41 is proven in Section [4.11 □ 
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3 Hyperbolic multifractal analysis 



We are now interested in the simultaneous analysis of local regularity properties and of 
anisotropic features of a function. To that end, we construct a new multifractal analysis, 
referred to as the hyperbolic multifractal analysis. Recall that in the isotropic case, the pur- 
pose of multifractal analysis is to provide information on the the pointwise singularities of 
functions. Multifractal functions are usually such that their local regularity strongly vary 
from point to point, so that it is not possible to estimate the regularity index (defined below) 
of a function at a given point. Instead, the relevant information consists of the "sizes" of the 
sets of points where the regularity index takes the same value. This "size" is mathematically 
formalized as the Hausdorff dimension. The function that associates the dimension of the set 
of points sharing the same regularity index with this index is referred to as the spectrum of 
singularities (or multifractal spectrum). The goal of a multifractal formalism is to provide 
a method that allows to measure the spectrum of singularities from quantities that can ac- 
tually be computed on real-world data. We extend this approach to the anisotropic setting. 
Let us first recall that, in the case where the anisotropy of the analyzing space is fixed, this 
has already been achieved: See [?] for anisotropic pointwise regularity analysis using Triebel 
bases and [?] for the corresponding anisotropic multifractal formalism. Here, we generalize 
these two previous works, providing a multifractal analysis which does not rely on the a priori 
knowledge of the regularity and takes into account all possible anisotropics. Note that for a 
fixed anisotropy, both formalims coincide: Indeed they are derived from wavelet characteriza- 
tions of the same functional spaces. Nevertheless, the formalism based on hyperbolic wavelet 
allows to deal simultaneously with all possible anisotropics, thus providing more useful algo- 
rithms for analyzing real-world data. In addition, the use of hyperbolic wavelet bases offers 
the possibility to define and synthesize deterministic and stochastic mathematical objects 
with prescribed anisotropic behavior. 

In Section [3. 11 the different concepts related to pointwise regularity are first recalled. An 
hyperbolic wavelet criterion is then devised in Section [3.1.21 Our main result, Theorem 13. H 
is stated in Section 13.1.21 and proven in Section HI Hyperbolic wavelet analysis is defined 
in Section 13.2.21 and the validity of the proposed multifractal formalism is investigated in 
Theorem 13.21 

3.1 Anisotropic pointwise regularity and hyperbolic wavelet analysis 
3.1.1 Definitions 

Let us start by recalling the usual notion of pointwise regularity (cf. ^35j for a complete 
review) . 

Definition 3.1 Let f be in Lf^^(M?) and s > 0. The function f belongs to the space C^^^ |/3 (xq) 
if there exist some C > 0, 6 > and Px^ a polynomial with degree less than s such that 

if \x - xo\ < 5, \ f{x) - Pccoix)\ < C\x - xo\' ■ I log(|x - xo\)f , 

where \ -\ is the usual Euclidean norm ofM?. If /3 = 0, the space C^Qg|o(a^o) is simply denoted 
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Anisotropic pointwise regularity is further defined as follows. Let P denote a polynomial 
of the form: 

(/3i,/32)eN2 

and let a = {ai, 02) be an admissible anisotropy. The a-homogeneous degree of the polyno- 
mial P is defined as: 

daiP) = sup{ai/3i + 02/32, a/^^^^Sj 7^ 0}; . 
Finally, for any t = (ti,t2) £ 1^^, the a-homogeneous norm reads: 

\t\a = +1*2!^/"' . 

We can now define the spaces C|'^j^g|^(xo). 
Definition 3.2 Let f G ^^^(M^), 

a = (ai, 02) such that ai + 02 = 2, | • |a, s > and /3 G M. 
The function f belongs to C^^^^^fj{xo) if there exists some C > 0, 6 > and P^o polynomial 
with a-homogeneous degree less than s such that 

if \x - xoU < ^1 \ f{x) - Pxo{x)\ < C\x - xol^ • I log(|x - Xo\a)\'^ . 

If P = 0, the space C|*j^g|o(2;o) is simply denoted C*'"(xo). 

The anisotropic pointwise exponent of a locally bounded function / at xq can be thus be 
defined as: 

(4) /i/,„(xo) = sup{s, /gC^'"(xo)} . 

The reader is referred to I10I.I36I for more details about the material of this section. 



3.1.2 An hyperbolic wavelet criterion 

As in the usual anisotropic setting (see |35j), the anisotropic pointwise Holder regularity of 
a function is closely related to the decay rate of decay of its wavelet leaders. The usual 
definition of wavelet leaders needs to be tuned to to the hyperbolic setting: 
For any (ii, j2, ^i, ^2), let A(ji, j2, A;i, ^2) denote the hyperbolic dyadic cube: 

^c-s , \ / • • r 7 \ , /ci A:i + 1 A;2 A;2 + 1 . 

(5) A = XijuJ2,kuk2) = [— , —^[ , 

and let cx stand for Cj^ j2,A;i,fc2- "^^^ hyperbolic wavelet leaders dx, associated with the hyper- 
bolic cube A, can now be defined as: 



dx = sup \cx'\ ■ 

A'CA 



For any xq = (a, b) G M^, let 



'i^jiM^o) - [ ' 2^2 ' 2^2 ^' 

(where [•] denotes the integer part) and 

dhjii^o) = sup |ca'| . 

A'CSAj-^ J2 (xo) 

The hyperbolic wavelet leaders criterion for pointwise regularity can now be stated as: 
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Theorem 3.1 Let s > and a = (ai, 02) G (M!j_)^ such that ai + 02 = 2. 

1. Assume that f £ C^'"{xo). There exists some C > such that for any ji,j2 G NU {— 1} 
one has 

(6) |d,,,,,(xo)|<C2-'"^^^^'^)^ 

2. Conversely, assume that f is uniformly Holder and that |^ holds, then f S C^^^^^2{xq). 
Proofs are postponed to Section 4. 

3.2 Anisotropic multifractal analysis 
3.2.1 Two notions of dimension 

In multifractal analysis, two different notions of dimension are mainly used: the Hausdorff 
dimension and the packing dimension, whose definitions are now recalled. 

The Hausdorff dimension is defined through the Hausdorff measure (see for details). 
The best covering of a set ii^ C M'' with sets subordinated to a diameter e can be estimated 
as follows, 

00 00 
n^,{E) = inf{ 1^,1^ Ec\jEi, \Ei\ < e}. 

i=l i=l 

Clearly, Ti^ is an outer measure. The Hausdorff measure is defined as the (possibly infinite 
or vanishing) limit Tif as e goes to 0. 

The Hausdorff measure is decreasing as 6 goes to infinity. Moreover, T-L^{E) > implies 

{E) = 00 if (5' < (5. The following definition is thus meaningful. 

Definition 3.3 The Hausdorff dimension dim.H{E) of a set E C is defined as follows, 

dim//(S) = sup{(5 : H^E) = 00} . 

With this definition, dim//(0) = —00. 

The box dimension (or Minkowski dimension) is simpler to define and to use than the 
Hausdorff dimension. 

Definition 3.4 Let E C M'^. If e > 0, let Ni;{E) be the smallest number of sets of radius e 
required to cover E. The upper box dimension is 

dimB(E) = limsup — ^ — ^. 

£^0 - log e 

The lower box dimension is 

dim^(ii;)=liminfi^^^^. 

£^0 -loge 

If these two quantities are equal, they define the box dimension dimsiE) of E. 

A significant limitation of box dimension is that a set and its closure have the same dimension. 
The packing dimension (introduced by Tricot, see [50]) has better mathematical properties 
(e.g., the packing dimension of a countable union of sets is the supremum of their dimensions). 
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Definition 3.5 The packing dimension dimp(£') of a set E CM."^ is defined by 



dimp(^) = mf{sup{dimB (Ei)} : E C [J Ei}. 



The following inequality holds for any set E C M"^, 

diuiniE) < dimp(£;). 

We now define the hyperbolic spectrum of singularities of a locally bounded function using 
the Hausdorff dimension. 

Definition 3.6 Let f be a locally bounded function and a = (01,02) G (MJ^)^ such that 
oi + 02 = 2. The iso~anisotropic-Hdlder set are defined as 

Ef{H,a) = {x£ M^/l/,c,(x) = H} 

where the anisotropic pointwise Holder has been defined in 

The hyperbolic spectrum of singularities of f is then defined as: 

d : (M+ U {00}) X (0, 2) ^ M+ U {-00} {H, a) ^ dimH{Ef{H, (a, 2 - a))) . 
3.2.2 The hyperbolic wavelet leader multifractal formalism 

It is not always possible to compute theoretically the spectrum of singularities of a given func- 
tion. A multifractal formalism thus consists of a practical procedure that yields (a convex 
hull of) the function d, through the construction of structure functions and the use of the Leg- 
endre transform. In the classical case, these formalisms are variants of a seminal derivation, 
proposed by Parisi and Prisch in the context of the study of hydrodynamic turbulence |43j . 
The hyperbolic wavelet leader multifractal formalism described below aimed at extending the 
procedure to where both anisotropy and singularities are studied jointly. 

Let us define hyperbolic partition functions of a locally bounded function as: 

(7) S{j,p,a) = 2~^^ Y: E d^nMM' 

where Tj[a) has already been defined in Section [2.21 with Eq. ([T]). 

From the definition of an anisotropic scaling function (or scaling exponents) 

(8) LOf [p, a) = hm mf — - — — — , 

j-s>oo log 2 J 

let us further define the Legendre hyperbolic spectrum: 

(9) Cf{H,a)= infjHp-ujf{p,{a,2-a)) + 2} . 



Qualitatively, the Legendre hyperbolic spectrum and the hyperbolic spectrum of singular- 
ities df{H,a) are expected to coincide, while the theorem below actually provides an upper 
bound relationship. 



13 



Theorem 3.2 Let f a uniform Holder function. The following inequality holds 
(10) V(F, a) G (M;) X (0, 2), df{H, a) < Cf{H, a). 

Definition 3.7 Let f a uniform Holder function, H > and a G (0, 2). If Eq. I jOj) simplifies 
into an equality, i.e., 

y{H, a) G (M;) X (0, 2), d{H, a) = Cf{H, a), 
then f satisfies the hyperbolic multifractal formalism. 

From an applied perspective, Eqs. ([7]) , ([8]) and ([9|) constitute the core of the practical 
procedure enabling to compute the Legendre hyperbolic spectrum from the hyperbolic wavelet 
leaders computed on the data to be analyzed. Practical implementation show preliminary 
satisfactory results, notably, for isotropic function, it is clearly observed that the measured 
Cf{H, a) does not depend on a. 

4 Proofs 

4.1 Proof of Theorem [272] 

4.1.1 Hyperbolic Littlewood-Paley characterization of Bp^{^) 

Let 6*0 G 5(M, M+) be supported on [-2, 2] such that 6lo = 1 on [-1, 1]. For any j G N, let us 
define 

0,- = 0o(2^-)-^o(2^-'-) • 

such that X]j>o^i(') = 1 is a resolution of the unity. 

Observe that, for any j > 1, supp(6'j) C {2^'^ < \^\ < 2-'+i}. 

Remark 4.1 In the following, the function 9q can be chosen with an arbitrary compact sup- 
port. It does not change the main results even if technical details of proofs and lemmas have 
to be adapted. It allows to chose the Fourier transform of a Meyer scaling function for 6q. 

Definition 4.1 1. For any j, £ >0, and any ^ = (^1,^2) G set 

For any ji, j2 > 0, the function 4>j-ij2 belongs to S{M?) and is compactly supported 
on {2^1 < 161 < 2^1+1} X {2^2 < 1^2! < 2^2+1]. p^^ther En>oEj,>o^h,n = 1 ^^id 
j2)(ii ,i2)eN2 is called an hyperbolic resolution of the unity. 

2. For f G 5'(M2) and ji, is > set 

The sequence ((Aj^ j2/)ii j2>o) called an hyperbolic Littlewood-Paley analysis of f . 



14 



In the remainder of the section, we are given a = (ai,Q!2) a fixed anisotropy satisfying 
ai + "2 = 2 and {f'j)j>o an anisotropic resolution of the unity. One then defines the following 
functions for any j > 0, 

(11) 9j= E 'I'n,.^ 

ji,j2e^j (a) 

where the sets Tj(a) have been defined in ([1]). 

Remark 4.2 Hyperbolic Littlewood-Paley analysis is used in the definition of spaces of mixed 
smoothness. We refer to J^T} / for a study of these spaces and to 14^1 for their link with tensor 
products of Besov spaces and their hyperbolic wavelet characterizations. 

We now provide the reader with the following hyperbolic Littlewood-Payley characterization 
of anisotropic Besov spaces: 

Theorem 4.1 Let s G M and {p,q) G (0, +oo]^. 

1. (a) If q < oo and 

(12) 

\ 1/9 

^^•qmax(l/p-l,0)+max(g-l,0) .^-/3g2i^9 ^ ||Aj,j2(/)|M < +00 , 

^i>o (ii,i2)er,(a) / 

thenfeB;2^^^^^,{R^). 
(b) If 

(13) max (i/p-i.o)+i . ^-P^^s ||A,,,,(/)||^) < +oo , 

2. (a) Ifq<^andf^ B'^;2^^^^^,{M?) then 

\ 1/9 

^^-1.^-/392.^9 ||A,,,,(/)||^ <+00. 

i>o (iij2)erj(a) / 

(h) then 



max 1 
i>o ' 



-^2^^ ||A,,,,(/)||J <+oo. 

{jij2)erj(a) 



The proof of Theorem 14.11 consists of several steps, beginning with 
Lemma 4.1 1. For any j >0 and any (ji, J2) '^j{ct)j one has 
(14) supp(^j") n supp((/)j, = . 
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2. For any j > and any £ ^ {j — + 1}, one has 

(15) supp(5rj") n supp((/??) = 



Proof. Point [J of the lemma is first proved, that is if {ii > Lmlx + 1 and £2 ^ -^max + 1) or 
ih < L^IL - 1 and £2 < 45n - 1). then supp(99") n supp(,/.,„,,) = 0. 



Indeed, if ^ G supp((^p, then ^ G Rf^^ \ Rf. Hence, if ii > L^rtL + 1 and £2 > L\^L + 1, 
one has for i = 1,2, |2~^*^j| < 2"^^^^^^^^^ < 2"*~^, by assumptions on £i,£2- Since ai or 
02 necessarily belongs to (0,1), one has ^ supp((/>£j^£2)- Hence, (pi^/^iO = 0- The same 



(2) 



approach leads to (C) = if ^1 < -^^min — 1 and £2 < — 1 and Point ([T]) of Lemma 14.21 
is obtained. 

Point ([2]) of Lemma 14.21 can be obtained similarly. □ 

From Lemma 14.21 an intermediate hyperbolic Littelwood Paley characterization of anisotropic 
Besov spaces is obtained. 

Proposition 4.2 Let {p,q) £ (0, +00]^, s,/3 G M. the two following assertions are equivalent: 
1. /G B^'",, ,«(M2). 



-(2) 



2. 



( 



(16) 



j>0 



{iij2)erj(a) 



1/9 



< +00 . 



Proof of Proposition \4-2\ 
Let us first show the 
this end, Point [H of Lemma 14.21 is used to deduce that for any j 



Let us first show that / G j^g|^(M^) implies Inequality (116p of Proposition 14.21 To 



,(iiJ2)eN2 



where 5" is defined by Equation (llip . Observe now that replacing the usual dilation with an 
anisotropic one gives an anisotropic version of Equation (13) in Section 1.5.2 in |51) . More 
precisely assume that we are given p G (0,+oo], 6 > and M G 5(M2). There exists some 
C > not depending on b nor M such that for any h G LP{M?) such that supp(/i) C {.^ G 
M^, supj l^il < 6"*}, one has 



(17) 



T-HMTh)U.m2) <C\\M{b^-)\\j,s 



LP 



where -ff| is the usual Bessel potential space and s > 2(1/ min(p, 1) — 1/2). 

Set now b = 2^ , M = ip'j and h = f . Since (p'j{2^°'-) = tpf, there exists some C > not 
depending on j such that for any p G (0, +00] and any / G LP(M^) 

< C\\ Yl ^Juhfhr = C\\{T-'gf) * fhr . 

(iij2)erj(a) 
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Then 
(18) 



1/9 



1/9 



B" 



,|log|/3 



I LP 



<c\Y,r^"in\{:F-^9'!)*f\\ 



LP 



which gives Point [1] of Proposition 14.21 

Let us now prove that if Equation (fTUj) of Proposition 14.21 holds then / G \ \og\0'' 
Point [2] of Lemma 14.21 gives for any j > 

9j f = 9j [^j-i + + ^j+i) f ■ 

Hence, Inequahty (fT7|) apphed with b = 2^ , M = and h = ('/'"-i + V'j + f'j+i)f gives the 
existence of some C > not depending on j nor / such that for any p G (0, +oo] 

II i^-'gf) * /IIl. < c||5,"(2^°•)llH|||(-F-V^l +-^-Vi +-^-V"+i) * /IIl. 

Since || • ||lp is either a norm or a quasi-norm (according to the value of p), there exists some 
C > such that 

II(-^-'5,")*/IIl. < c |b,"(2^"-)lki 

(IK-^'Vi-i) * fh. + IK-^-V") * fh. + IK-^-Vi+i) * /IIl.) . 

Let us first bound ||(7°(2''"-)||//|. To this end, observe that 

^bi(2^"-)](6 = 2"^(°^+°^)g;(2-n) = 2-2^-5;(2-n) 

^ 2-''^+^2-2j?i(2^i--'°i^i)^i(2-''^-"2i^2) . 
0ij2)erj(Q:) 

Hence 

2 



Il5^(2^"-)ll?/^ 



R2 



(1 + iei 



2ns 



.(iij2)erj(Q:) 



de 



< 



/ (i + iepr 



.(iij2)erj(a) 

Since 9i £ 5(M), for any > 1 there exists some C > such that 

Cm 



de 



l^i(C)l < 



(i + ICI)2^ ■ 



Finally 



||5,"(2^"-)ll^| < Cm I {l + m' 



< Cm 



[ (1 + ICI 

Jr2 



E 

(jij2)erj(a) 

2\s 



2ii+i2-2j 



(1 + |2Ji-i"i ^i|)2Af . (1 + |2i2-ja2 ^2\)2M 



d^ 



(iij2)Grj(a) 



(2i"i-ii + |^i|)2A^ . (2i«2-i2 + |^2|)^^ 



de. 
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By the inequality 

(a + 6)^ > a max(6, 1) , 

valid for any a > 1, 6 > and applied successively with a = 2^°^^"^^ and b = a = 2^°'^~^^ 
and b = \^2\, it comes 

||<(2^"-)llk^ < Cm [ + 



. , max(l,|eiirmax(l,|6l^*^ 

With a M sufficiently large it follows that 



dC 



sup(||5^(2^"-)k|) <+oo 



Going back to an upper bound of || ( J-" I * fWip, there exists some C > such that 



and 



1/9 



(19) |5]r^'2^'"'ll(-7^-'5,")*/ll 



q 

LP 



<c^II/IIb- , 

p,9,l log r 



9 

LP 



the last shows that if ([T6]) holds then / G 5j^n^g|^(M2) 



Proof of Theorem \4- 1\ Let us first recall that: 

• For any p £ (0, +oo], n G N, and (/i, • • • , /„) G LP(R2)" 

(20) ll/i + • • • + fnWLp < n^-^i/^'-i'O) (ll/ill + . . . + ll/JI) 

• For any g G (0,+oo), n G N, and (ai,--- ,a„) G (1^+)" 

(21) (ai + • • • + a^y < (af + • • • + a^) . 

Let us now prove the first point of the theorem in the case where q ^ oo. For this, assume 
that ([12]) holds and let us prove that / G 5^'^n^g|^(R2). By Inequalities ([20]), ^ and the 
fact that Card(rj(a)) < Cj there exists C > such that 



Hence, 



(jij2)erj(a) 



J>0 



LP 



(iij2)erj{a) 



9 

LP 



(iij2)erj(Q) 



9 \ 1/9 



LP/ 



1/9 



vi>o (iij2)er,(a) 



9 

LP 
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It proves that if (jl2p holds, one has 



1/9 



9 

LP 



< OO . 



Finally, by Point (1) of Proposition 14.21 it comes that f €z 
We now deal with the case q = oo. In this case, we have 



,9,1 log 1/3' 



max j-'^2^' 



(iij2)erj(a) 



LP 



< Cmaxj--(i/P-i,o)^-/32i'' ^ ||A,„,J|Up . 

(iij2)6r,{a) 



Hence if m holds, / E i?;;;^u^g|,(R2). 

To prove the converse assertion, let us assume / G Observe that for any 

j > and any (ji, j2) G rj(a), one has 



(Pji,j2f = <Pji,j2 {9j-i + 9j + 9j+i) f ■ 

Remark that ^j^jjl^"'"") is bounded in i7|(R^) independently of (ji, ^'2) S rj(a). Hence, by 
(fT7|) . there exists C > not depending on j nor / such that for any (ji, 72) S ^jic^) 

Again, two cases have to be distinguished according whether g 7^ 00 or g = 00. 

Let us consider the case q < 00. Observing that Card(rj(a)) < Cj, we deduce that 



E ii(-^"Vi„i2) * Wlp < Cj E wi^-'ar) * /II 

(iiJ2)Grj(a) l=j-l 



9 

LP 



So 

(22) E E ii(-^"V,„,2) * /Hi. < E • i-'j-''^2^-^^ii(^-i5,") * /II 

i {jij2)erj{a) i 

Since in addition the function / is assumed to belong to one has 

Y^j-^'^2n\{^-'9j) * Wlp = E^' • r'r^'2n\{^-'9j) * /lli. < oo , 



9 

LP • 



which directly yields the required inequality using 
In the case q = 00, we have 

max ||(^-Vji,,.)*/|lip < C max ||(^-V)*/IIl- 
0ij2)erj(a) ^=j-i,jj+i 

which leads for some C > to 



max 1 J -^2^^ max || (-F" * /||l. ) < Cmax (j-''2^m{T-'g-) * /Up) , 
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that is 
(23) 

max (max(^, ^)) 2'""^(^'^^1(^- V,,,,.) * fhr < Cms.xj-^2^'\\{T-'gf) * fU, . 
Since in addition / is assumed to belong to -B*'^ it comes 

ma^j-^2^^\\{T-'gf)*f\\LP<oo. 

Finally, the required conclusion is obtained by an approach similar to the one used for the 
previous case. 

4.1.2 The special case p = 2 

In that case, an exact hyperbolic Littlewood-Paley characterization of anisotropic Besov 
spaces is provided: 

Proposition 4.3 Let f G 5'(M^), s > and q G (0,+oo]. The two following assertions are 
equivalent 

q q 

i>o V{jij2)er,(a) J j>o \Uuh)er,(a) J 

+00. 

In particular the following exact hyperbolic Littlewood-Paley characterization of anisotropic 
Sobolev spaces can be stated: 

Theorem 4.4 Let f G 5'(M^), s > and q G (0,+oo]. The two following assertions are 
equivalent 

(n) E r^P2^^' E Un,J\\h = E E W^nADWh < 

To prove Proposition 14.41 l^t us first precise the relation between anisotropic and hyperbolic 
resolutions of the unity. 

Lemma 4.2 For any j > 0, the following inequality holds on 

(24) ^'j<9j= E <l^n,n<Vj-i + ^'! + ^'j+i- 

(jij2)er,-(«) 

Proof. Let us first observe that 

0ij2)eN2 
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To get the left hand side of inequahty (p^ . we just have to prove that if ^ = (^1,^2) £ 
supp((^"), one has ^jijalO = if (ji,i2) ^j, that is if {h > L^L + 1 and £2 > imix + 1) 



or (^1 < - 1 and ^2 < L^l - !)• 



(2) 



Indeed, if ^ G supp((^j"), then ^ G \ i?^". Hence, if li > L^^L + 1 and h > ^mL + 1, 
one has for i = 1,2, |2~^»^j| < 2"^^^^^^^^^ < 2"*~^, by assumptions on £i,i2- Since ai or 
02 necessarily belongs to (0,1), one has ^ supp((/)£j^£2)- Hence, </'£i/2(0 ~ 0. The same 
approach leads to 0^i,^2(C) = if ^1 < — 1 and ^2 < -^mL ~ 1- The left hand side of 



(2) 



inequality is obtained. 

Let us now prove that the right hand side of inequality (I24p holds. It comes from the 
obvious equality 

(jij2)er, i>o 

and if ^ e supp((/>ji for some (ji,i2) e rj(a) then for any I e {j - + 1}, cpK^) = 0. 

□ 

Before proving Proposition 14.41 l^t us give some background about quasi-orthogonal systems. 



Definition 4.2 Let a Hilbert space H and (•,•) the associated scalar product. A system 
{fk, k £ Z} of H is said to be quasi- orthogonal if there exists some £ £ N such that 



(25) 



y{k,k')€z\ {\k'-k\>e^ {fk,fk') = o) 



Lemma 4.3 Let H be a Hilbert space and \\ ■ ||,(-, •) the associated norm and scalar product. 
Let {fm, m £ 7j} a quasi-orthogonal system of H and let £ G'N satisfying \25\) . Then 



(26) 



<{2£+l)Y,\\fm\ 



Proof. Observe that for any m! G Z, {fm, fm') = except ii m' — £ < m < m' + £. Hence 



E -^"^ 

mall 



< E T.\(fmJm')\ 



< 



m'+e 



< E E 



m'fc£, m=m' 



m'+e 



< E^/2ZTT E II/" 



\m=m' —£ 



^ m I 



^^2^ E E II/" 

Vm'eZ ^ 

< (2^ + 1) Yl II/" 



Km'alj \m=m' 
|2 



E II/™' 

Vm'eZ 



m'GZ 
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□ 



Proof of Proposition 14.41 

By Plancherel Theorem and by Lemma |4.2[ one has 



and 

/ 9]{0\m\'dC < / M-i)' + (^i )' + (^i )'] 1/(01' < ||A^_i/||i. + ||A°/||i. + ||AJ_,i/||i. , 

where gj is defined by (j24|) . Proposition 14.41 is then a straightforward consequence of the 
following lemma : 

Lemma 4.4 There exists some C > such that for any j > 0, one has 

(27) Un,nf\\h<\\9jf\\h<C Un,jJ\\h- 

0ij2)erj(Q:) (ii,i2)erj(a) 

Proof. Since supp((^jjj2) n supp((^m,i,m2) = if max(|mi — ji|, |m2 — ^2!) > 3, the system 
(<^iij2/) is quasi-orthogonal. Hence, by Lemma there exists some K > such that 



0i,i2)er,(a) 



L2 



< ^ E ii'^.i-.2/iii2 

(iij2)er,(a) 



For the converse inequality, observe that each term {4>ji,j2f^ J^^^ ~ 'f'jiJiiO'Pj'-^j'^iOlfiOl'^d^ 
is positive. Hence 

Y ll'^ilJ2/lli2 < Y i^juj2f,<l>jld2f)L2+ Y {<Ph,hfAj[,j'J)L^ 

{h,j2)&j{a) (iij2)er,(a) (iij2)^(i(j^)erj(Q) 

2 



Y '^il'J2/ 

(ji,i2)erj(a) 



L2 



□ 



4.1.3 Proof of the hyperbolic w^avelet characterization of anisotropic Besov spaces 

Let us first consider the general case where {p,q) G (0, +oo]^,/3, s G M and a = (01,02) a 

^) are defined as the collection of functions 



fixed anisotropy. Intermediate spaces £ ' 



fofS 



such as 



A norm on £^'",, ,«( 
p,q,\ log 1'^ 



P>'?>|log|'' 

Yr^"^''' E iiA,„,,/ii^ < +00. 

i>o (jij2)er,(Q) 
is defined as follows 



Piij,! log 1*^ 



0ij2)er,(a) 



such that the embeddings 
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• if g < oo 



^p,'(?,|log|/3-m''''(l/p-1.0)-max(l-l/q.O) (I^ ) ^ p,q ,\ log {1^ ) "^p.g, | log | -S + V? ) ' 

• if g = oo 

'^p,'oo,|log|/3-m=->=(l/p-l.0)-l(^^) ^ -^piooJlogl/S^^^) ^ '^p.g.llogl/'^^^) • 

are an exact rewriting of Theorem 14. 1[ 

In the special case where p = 2, we proved in Proposition 14.41 that H^^^^^p{M?) = 

^2^1 log \^ (^') |, m coincide. 

In the following proposition, an hyperbolic wavelet characterization of spaces £plq{M?) 
is given. Combining Proposition 14.51 Theorems 14.11 and ?? directly implies Theorems 12.21 
andlM 

Proposition 4.5 Let {p,q) € (0,+oo]^, s,fi ^ M^. The following assertions are equivalent 

Let us prove Proposition 14.51 The equivalence between assertions (2) and (3) holds since for 
any (ji, ja) G r,(a), one has max(l^, ^) + 2 - 2 < j < max(^, ^) + 2 and UT,- = Ng. The 
crucial point is the equivalence between assertions (1) and (2) . 

Proof of implication (1) =^ (2) of Proposition 1475) 

The proof of this implication relies on the following sampling lemma which is a adaptation of 
Lemma 2.4 of [27J in the case of rectangular support. 

Lemma 4.5 Let p G (0,+oo] and j = (ji,j2) G Nq. Suppose g G and supp(^) C 

161 < 2^1+^ and l^l < 2^2+1}. Then there exists C >0 such that 

\(fei,fe2)ez;2 V / y 

Proof. Let ip £ ^(M^) be such that supp(V^) C {6 max(|6|, |6|) < vr} and = 1 on [-2,2]^. 
Set ipj{x) = 2J'i+J2^(2^ixi,2^22;2). One has ^ = 1 on [-2^^+'^ ,2^^+^ x [-2^2+1,2^2+1]. 

By assumption supp(5) C [-2^1+^2^1+1] x [-2^2+1^2^2+1]^ that for any x = {xi,X2) G 

and any fixed y = (yi, 1/2) G 

g{x + y) = {i;j*g){x + y) = (2vr)-2 / / V5(05(Oe'"«e'^«de . 

y5i=-2Ji+i jg2=-2J2+i 
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Denote hj the periodic extension of i/^j with period 2^^^^7r for each variable {i = 1,2). 
One has 



(28) <7(^ + y) = (2vr)-2 / / (/i,(e)e^"« ) fe^e^^^ ) d^ 



Using an expansion of hje^^^ in two dimensional Fourier series, it comes 

/r,(e)e-« 



(^i/2)ez- 



2^1 + ^77 /.2-'2 + l7r 



= 2-in+j2) V'^.(x-2-^'^)e'2~^^^i«iei2-^2£26 ^ 

(£i,fe)ez2 

where for j = (ji,i2) and £ = (^1,^2), the notation 2~H = {2~^^li,2~^'^(.2) is used. Replacing 
hj{^)e^^^ with the last sum in Equation (f28l) yields that for any x = (xi,X2) S and any 
fixed y = (2/1,^2) e 



5(2; + y) 
2-(ii+i2) 



47r2 



y / / V'i(a;-2-^'^)e^2-n^,aei2-^^26(^(^)eis/C)de 



Hence for ah y G Xj.j^MM = [2-^'ki,2-^^ {h + 1)) x [2-^2/02, 2-^2 (/c2 + 1)) 
sup \g{z)\ < sup \g{x + y)\ 

\zi-2-hki\<2-h,\z2-2-hk2\<2-^2 \xi_\<2-h ^\x2\<2-^2y/2 

< 2-(^'i+-'2) \g{2-H + y)\- sup IV'i(a;-2" 

(^i,^2)eZ2 max(2W|xi|,2J2|a:2|)<x/2 

where the last inequality follows from the fast decay of tp. Take M sufficiently large and use 
either the p triangular inequality either the Holder inequality according whether p G (0, 1) or 
p G [l,+oo]. Hence, one has 

|5(2-^'^A;i,2-^2^2)|P< |g(z)|f < C72-(^-i+^2) ^ |g(2-^£+y)|P- ] 

\zi-2'nki\<2-n,\z2-2~^2k2\<2-n (4/2)6X2 U + I^U 

for some M' > 1. An integration over y G Xji,j2M,k2 leads to 

2-0i+.2)|^(2-^i^^^2-^2^^)|P< Y JYTJmW [ \9iy)\'dy 
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and a sum over k ^ 7? gives 



k {fi/2)ez2 



l5(y)rdy 



which ends the proof of Lemma 14.51 

Now, observe that cj-^^j^^kiM ~ ^iij2/(2~"'^^i5 ^"-'^ ^2). By Lemma [4.51 applied to g 
Aj.jJ E 5(M2), one has ' ' 

(A;i,A;2)eZ2 (fci,fe2)GZ2 

which is the desired wavelet characterization. 



□ 



Proof of implication (2) =^ (1) of Proposition 1475] 

To obtain the converse implication, the same approach as in the proof of Theorem 3.1 of [27j 
is followed. 

Since 4>jij2 '^mi,m2MM both defined as a tensorial product, Lemma 3.3 of [27] 
can be applied: there exists some C > such that for any a > and for all x = (xi, X2) G 
one has 
(29) 

2-(|ji-mi| + |i2-m,2|){Af+3) 
l'/'ji,J2*V'mi,m2,fci,fci(2;)| < C 2^ni(n^m^)\xi - 2-'"l A^i | )" (1 + '2}r,i(n.m2)\x2 - 2-^^k2\Y ' 

where M denotes the number of vanishing moments of the wavelets. 
A lemma analogous to Lemma 3.4 of [27] is now proved: 

Lemma 4.6 Let p G [1, +oo], ^i, ^2) m2 integers such that ii < nii and I2 < m2. We are 
also given some functions gki,k2 satisfying the following inequality for some C > : 



(30) Vx = (xi,X2)GMM5fc,,fc,(x)| < 



C 



(1 + 2^i|xi - 2-™i/ci|)2(l + 2^2|x2 - 2-™2/c2|)2 ■ 



Set 



Then 



(31) 



F = E dkiMdkiM 



,fc=(fci,fe2)ez2 



i/p 



Proof. By definition of the L^-norm, one has : 



l-f^ Wlp 



E dk^MdkiMix) 
fc={fci,fc2)ez2 



dx 



hi — (hi hi \rz-7/% •' ^ 



k' = {k'^,k'2)ei? " ^mi,m2,*i,fe^ 



E dk^MdkiMi^ 

fc=(fci,A:2)eZ2 



dx 
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where the hyperbohc dyadic cube Xmi,m2,k[,k'2 defined in ([5]). Observe now that, by the 
usual triangular inequality and by inequality (j30p . there exists some C > such that for any 
{ki,k2) G Z2, {k[,k'2) E Z2 



sup 

m-^ ,m2 , fe^ 



fc=(fcl,A;2)eZ2 



^ ^ I^fci,fc2 



fe=(fci,fc2)ez2 ii«-i,2v ' I 2 



Hence one has 



(fci,fc^,)ez2 \(fci,fc2)GZ2 V ' 11 ii7 V ^ 12 ^ly 

Let us recall the usual convolution inequality, valid for any sequences s, s' in (PiW?) for p>l, 

II 'l|P ^ II \\P II 'I|P 

Applied to s = \dk,M\ and s' = (1 + 2^i-"i | fc; - ki\)-'^{l + 2^2-™2|^^ _ k2\y^, it gives 



1 



(l + 2^i-™i|A;;|)2(l + 2^2-™2|A;^ 
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^2 



^(fci,fc2)ez2 y \(fci,fc^,)6Z2 

Recall now the classical result : 

\^ \ <- (<-<nmi-£i 2^2-^2 

.'=(.t^)eZ2 (l + 2^--|^'ll)2(l + 2^2-rn2|fe.|)2 - 

Hence 

||F||^p < C72-{'»i+™2)2("^i-^i)P2(™2-^2)p I ^ |4i,fc2r 

\A:=(fei,fe2)eZ2 
P 

y ^ ' 

....,.t^,a.n„,.(i + 2'.-".|t:i)^ 

which directly yields the required result. It ends the proof of Lemma |4.6[ □ 

Let us now go back to Implication (2) =^ (1) of Proposition I4.5[ Two cases are considered: 

p £ (0, 1) and p £ [l,+oo]. 

Let us first assume that p £ (0,1). 

We have to bound || A^^ Jj/IIlp = Il'/'jij2 * IWlp- Observe that 

4^31^2 * / ~ ^ ^ ^ ^ Cmi,m2,ki,k2 i.4'jl,j2 * ^r?ii,m2,fci,A;2) 
Til ,7712 k-i,k2 

By the p-triangular inequality, it comes 

Vx= (X1,X2) G |</'jij2 ^ X] X] |Cmi,m2,fci,fc2ri(</'jlj2*V'mi,m2,fci,fc2)(2;)r 

mi,ni2 fci,fc2 
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By Inequality (p9]) . for all x = {xi,X2) G M^, one has 



mi,rn2 k\,k2 



2^P(lii^'"i|+b2-m2|)(A/+3) 



(1 + 2i°f(jl'™l)|xi - 2"™l/ci|)P°(l + 2i"f02'™2)|a;2 - 2-™2fc2|)P" 



An integration over implies that : 



Hence 



mi,m2 ki,k2 



-I3q 



■p(lii-mi| + li2-m.2|){Af+3) 



"p.gJlogl/S \ Oil «2 

J1J2 



LP 



< 



EE 



-nii,m2,-r ll£p 



ji,i2 \mi,m2 

Set for any t G M, = max(t, 0) and 



sgn(t) = < 



' 1 if t > 0, 
if t = 0, 
-1 if t < 0. 



Observe now that for any integers j, m 

m — (m — < J < (j — + m 
and that for any integers Ji,j2)"ii,m2 



max(^, ^) 



1 



, 1 / , — < maxf — , — ) < max( , ) 

max{l^il,l^) 



1 + max 



ai 



ai 



(except in the case mi = m2 = which can be treated separately). Hence 

/ \ g/p 



< 



p. 9,1 log 1*^ 



EE 

J1J2 \mi,m2 



with 



and 



/ .mi m2. 

max( — , — 
V ai a2 



/ mi mo , 

„psmax( — ^, — ^)|| IIP 
^ ± ^ <'mi,m,2,-r ' 



^Jlj2 - 2 



-P{|il| + |J2|)(M+3)J-^ _^ max(^^i^, (-^2)+ ^j-/3poSgn(s)psmax(% 



(J2) + 



If g/p > 1 Young's inequality can be applied, which states that for any sequences s, s' , 



* A\£1/P(I?) < ll'S|l^9/P(Z2)ll'5'll£l(Z2) , 
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whereas if q/p < 1 the usual (g/p)-triangle inequahty and the usual inequality ||s * s'||£i(22) < 
||'S'||£i(z2) valid for any sequence s, s' are applied. In any case, the following inequality 
is obtained 



9 

p.g,| log \^ 



f sr^ f , mi 7722 A ^^„g.max(^,^)|| 

< > max( — , — 2^ "1 "2 c„i,„2, 



J2 2-^'(^-^+^-2)(^+3) (max(^, ^) j 2^'^"''^^(^'^) 



X 



If the wavelets have sufficiently vanishing moments, we get that 

-I3p 



P,<1,\ log 1^ 



which is the required result. 

We now consider the case p S [1, +oo]. In this case, observe that 

ki,k2 

with 

yki,k2 — \r'ji,j2 fmi,m2,ki,k2) ) 

and 

J , _ 2"(lii-™ll + li2-ni2|){Af+3) . 

"'Kl,K2 "-Jl,J2,Kl,«:2 • 

We set £1 = inf(ji,mi) and £2 = inf(j2,m2). Lemma [461 gives 

HA flU ^ /^9-p{|ii-"ii|+li2-m2|)(M+3)9-(mi+m2)/p||„ ip r,mi-£i9m2-^2 

II JI1J2J ll-tj*' — ^ ^ \\^'rni,m2,-,-\tP^ ^ 

Again two cases q < 1 and q > 1 are distinguished and the same approach than in the case 
p G (0, 1) is followed. It leads to the required conclusion. 

4.2 Proof of Theorem SU 

First a two-microlocal criterion is proved. 

Proposition 4.6 1. Assume that f G C^'°'{xq). Then there exists some C > such that 
for any {juj2MM) G (NU{-1})2 x Z^, 



— a 

2J1 



32^ 
,2 "2 + 



2J2 



(32) |cjij2,fei,fc2l < C'min(2 + 

2. Conversely, assume that f is uniformly Holder and that i32\) holds, then f G C^^^^p{xo). 
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Proof. Let us first assume that / G C^'°'{xo) with xq = {a,b). Assume that ji ^ —1 and 
j2 7^ —1. By definition of the hyperbolic wavelet coefficients one has 



Since '0 admits at least one vanishing moment, the two following equalities hold 

(33) Cj,,,,,k,M = I (/(^i,^2) - P.„(a,X2))^(2^'^xi - ki)^{2^^X2 - k2)dxidx2 
and 

(34) Cj,,j,,k,M = 2''^'' / {f{xuX2) - P,„(xi,6))^^(2^ixi - ki)^{2^^X2 - A;2)dxidx2 
Equality ([33]) and the assumption / G ^'"(xo) imply that 

\cn,nMM\ < 2^'^'' / W " <\^{'^''^i " fci)^(2^'^X2 - A:2)|dxidx2 



< 2^'i+^2 



Xl - 

2-'i 



+ 



2^1 



V'(2^i2;i - A:i)V'(2J2x2 - A:2)|dxidx2 



We now set ui = 2^^xi — ki, U2 = 2^^X2 — k2 and deduce that 



Jl£ 



Cji,j2,fci,fc2| < I 2 "1 / |ui|^/°i|?/'(ui)V'(u2)|duidu2 + 



2^1 



s/ai 



|V'(ui)V'('"2)|dtildU2 



Hence for some C depending only on ^, s and a one has 

ki 



{cjujiMMl ^ <^(2 "1 + 
A similar approach yields that 

K,j2MM \ < C{2'^ + 



2-?i 



2^2 



s/ai 



s/a2 



This shows that (I32p can be read as a necessary condition for pointwise regularity of function 
/• 

Let us now prove the converse result. Assuming that (I32p holds, the aim first consists in 
defining a polynomial approximation of / at xq. To that end, a Taylor expansion is used to 
investigate the differentiability of / at xq. Let us define fj as: 

(ii,i2)er,(a) (fci,fc2)ez2 
where the notations are the same as in the proof of Proposition 12.31 One has 



i/.(^)i < E E 

(iij2)Grj (fci,fc2)ez2 



min(2-J'i^/"i + \^- a|^/"i,2--'2^/°2 + |^ _ a| 



(l + |2iixi-A;i|)^(l + |2i2x2-A:2| 



< 



2"-''" + l|j - X2\'^"' + \X2 - 



(l + |2iixi-A:i|)^(l + |2i2x2-fe2|)^ 

J1<J fci,fc2 J2<J 



2-^^ + 1231-^1 



1 + Ixi — a 



s/«2 



+ |2Jixi -A;i|)^(l + |2J2x2- A:2 
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Then 

(35) \Mx)\<C{j2-^'+j\xi-a\'/''^+j\x2-b\'/'^^) . 
In the same way, if /3 = (/3i, /32), one has 

(ii,i2)er, (fci,fc2)ez2 v i v i ^i; 

Then 

(36) \d^fj{x)\ < C2J'(/3i"i+/32a2)(2-i^ + _ a|«M + |a;2 - ^ 

So, the function / is /3-differentiable at xq provided that /3iai + /32a2 < s. The Taylor 
polynomial of / at xq is defined by 



/9iai+/3202<s 

and 

Pxq{x) = ^^-fj>o(^) ■ 

We shall now bound \f{x) — Pxg{x)\ in the neighborhood of xq. Recall that / is assumed to 
be uniformly Holder, namely there exists some Eq > such that / G C^o(M?). The inclusions 
between Holder spaces with different anisotropies (see [53|) leads to the existence of Eq such 
that / G C^0'"(R2). Set Ji = [aJ/eo]. Observe that 

Ji 

i/(x) -p,„(x)i < ei/,(x) -p,-,„(x)i + \fji^)\ + E + Ei^-^-'-o(a;)i ■ 
j<j j=j+i j>Ji j>j 

Let us now bound each term of the right hand side of this inequality. 

We first deal with the term corresponding to j < J. In this case we shall use an anisotropic 
version of Taylor inequality which can be found in |19j . |26j and recalled in [9]. It gives the 
existence of some C > such that 

- Pj,xo(.x)\ <C Yl 1^ - ^012^^^+"^''^ _ sup \d^f,\ . 

/3i+/32<fe+l,ai/9i+«2/32>s z-{zuZ2)m^ 

with k = [max(s/ai, s/a2)]. The bound ([36]) implies that there exists some C > such that 

\fj{x)-Pj,x,{x)\ <C J2 |x-xo|2i'^^+"^'^^2^'('^i"i+''2a2)(2-i^+|xi-a|^/"i+|x2-6|'/"^) 

/3i+/32<fc+l, ai/3i+02/32>s 

Hence, 

^|/j(x)-P,>o(x)| <C Y |x-xo|:^^^i+"^'^2(2-^(^^"^+''2"2-^)+2-^(^i"i+''2a2)|^_2.Q|s) 

j<J I3i+l32<k+l,ail3i+a2l32>s 
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Since \x — xqIo < 2 it comes 

(37) <C7|x-xo|^ 

j<J 

Let us now bound the sum Yl'j=j+i \fj{^)\- (133) and the definition of Ji which depends 
on J, one has 

Ji Ji 
i=J+i i=J 

To bound the sum X^j>j^ l/j(3^)l t^^^ uniform regularity of / is used, leading to 

(39) \fj{x)\ < 02--^^^" < C72-^" 

j>Ji 

the last equality following from the definition of Ji. 

Finally, by (p6]) . the sum Ylj>j \ Pj,xo{x)\ can be bounded. Indeed, for some C > 0, one 

has 



02 



j>J /3iai+/3202<s i>J /9iai+/32a!2<s i>J 



Since |xi — a| < |x — xqIq^ < 2 -^"^ and |x2 — 6| < |x — xqI"^ < 2 -^"^ 



comes 



(40) E I^J>o(^)l < <^ E 2--^('^i"i+'^2a2) ^ 2J(''i"i+^2a2-s) < (^2--^^ . 

j>J /3ioi+/92a2<s i>J 

Finally, Inequalities dST]), ([MD, dM]) and (gOD yield that / E C|''j'"g|2(2;o). □ 

Theorem 13.11 is a straightforward consequence of the two-microlocal criterion and of the 
following lemma: 

Lemma 4.7 The two following properties are equivalent: 
(i) Inequality holds, 
(ii) Inequality ^ holds. 

Proof. Assume that (p2]) holds. If A' C 2>\j-^^j^{xQ), then 

j'l > jiJ'2 > h , 

and 



Condition ([52]) implies 



- a| < 2.2"-'i and 1-^ - 61 < 2.2"^'2 
2-?i 2-?2 



|cA'|<min(2 °i , 2 °2 ) = 2 ^qi'q2^ 
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Set 



and 



Conversely, assume that ([6]) holds. Let A' = X{j[, k[, an hyperbolic dyadic cube. 

ii=sup{4,2-^i + A-a| <2-^i} 
2-^1 

J2 = sup{£2, 2-^U A-6|<2-^^} 

2-J2 



We have A' C 3Xj-^j^{xo). Since ([6|) holds one has 

— 111 _22£ _iLf 

|ca'| < min(2 °i ,2 °2 ) < (711^1(2 °i + 



2^-i 



,2 "2 + 



^2 
2^ 



that is [32] holds. 



□ 



4.3 Proof of Theorem [372] 

The proof of Theorem 13.21 is based on the two following lemmas, analogous to Propositions 7 
and 8 of [35]: 

Lemma 4.8 Set a = (a, 2 — a) and define 

Let p > and s G {0,uj{p, a)/p]. Then for any H > s — 2/p 

dimH{G{H, a)) <Hp- sp + 2 . 

IfH<s- 2/p, dimH{G{H, a)) = 0. 

Lemma 4.9 Set a = (a, 2 — a) and define 

B{H,a) = {x G / e C^'°(x)} . 

Let p < and s G (0, uj{p, a) /p] . Then 

dimniBiH, a)) < dimp{B{H, a)) < Hp- sp + 2 . 

The proof of Lemma [4.8l in the case H > s — 2/p is exactly the same as this of Proposition 7 
of [35], except that the set Gj^h are replaced with the sets 

Gij,H,a) = {X = Xiji,j2,h,h), (ji,i2) e r,(a),|dA| > 2-^'^P} . 

Lemma 14.81 in the case H < s — 2/p, comes from the hyperbolic wavelet characterization of 
anisotropic Besov spaces stated in Theorem 12.21 and the Sobolev embeddings which can be 
proved in the anisotropic case as in the isotropic one (see [53]). 

The proof of Lemma 14.91 is exactly the same as this of Proposition 8 of |35j , except that 
the set Bh are replaced with the sets B{H,a). 

Lemmas l4.8l and l4.9l then imply Theorem l3.2l since for any a = (ai, 02) such that 01+02 = 
2 one has 

EiH,a) C {nH'>HGiH',a)) D {Uh'<hB{H' ,a)) . 

Acknowledgements. We warmly thank Florent Autin and Jean Marc Freyermuth for 
many stimulating discussions about applications of non parametric statistics to the analysis 
of anisotropic textures as well as Laurent Duval for giving us very interesting additional 
references about hyperbolic wavelet analysis. 



32 



References 

[1] P. Abry, S. Jaffard, and W. Wendt, Irregularities and scaling in signal and Image 
processing: Multifractal analysis, "Benoit Mandelbrot: A Life in Many Dimensions" 
World Scientific M. Frame Ed., 2012. 

[2] H. AlMAR and I. Gomez, Parabolic Besov regularity for the heat equation., Constr. 
Approx., 36 (2012), pp. 145-159. 

[3] F. AuTiN, G. Claeskens, and J.-M. Freyermuth, Hyperbolic wavelet thresholding 
rules: the curse of dimensionality through the maxiset approach., Submitted, (2012). 

[4] A. Ayache, Hausdorff dimension of the graph of the Fractional Brownian Sheet., Rev. 
Mat. Iberoamericana, 20 (2004), pp. 395-412. 

[5] A. Ayache, S. Leger, and M. Pontier, Les ondelettes d la conquete du Drap Brown- 
ien Practionnaire., C. R. Math. Acad. Sci. Paris, 335 (2002), pp. 1063-1068. 

[6] A. Ayache, F. Roueff, and Y. Xiao, Local and asymptotic properties of Linear 
Fractional Stable Sheets., C. R. Math. Acad. Sci. Paris, 344 (2007), pp. 389-394. 

[7] , Linear Fractional Stable Sheets: Wavelet expansion and sample path properties.. 

Stochastic Processes and their Applications, 119 (2009), pp. 1168-1197. 

[8] A. Ayache and Y. Xiao, Asymptotic properties and Hausdorff dimensions of Fractional 
Brownian Sheets., J. Fourier Anal. Appl., 11 (2005), pp. 407-439. 

[9] H. Ben Braiek and M. Benslimane, Baire generic anisotropic multifractal formalism 
in criteria in anisotropic function spaces. Submitted, (2011). 

[10] , Directional and anisotropic regularity and irregularity criteria in triebel wavelet 

bases. Submitted, (2011). 

[11] M. Benslimane, Multifractal formalism and anisotropic self similar functions. Math. 
Proc. Camb. Philos. Soc, 124 (1998), pp. 329-363. 

[12] G. Beylkin, Wavelets and fast numerical algorithms, in Proceedings of Symposia of 
Applied Math., vol. 47, 1993. 

[13] G. Beylkin, R. Coifman, and V. Rokhlin, Fast wavelet transforms and numerical 
algorithms I, J. Comm. Pure Appl. Math., 44 (1991), pp. 141-183. 

[14] H. Bierme, M. Meerschaert, and H. Scheffler, Operator scaling stable random 
fields., Stoch. Proc. Appl., 117 (2009), pp. 312-332. 

[15] A. BONAMI AND A. Estrade, Anisotropic analysis of some Gaussian models., J. Fourier 
Analysis and Applications, 9 (2003), pp. 215-236. 

[16] M. BOWNIK, Anisotropic Hardy spaces and wavelets, Mem. AMS, 781 (2003). 

[17] , Atomic and molecular decomposition of anisotropic Besov spaces. Math. Z., 250 

(2005), pp. 539-571. 



33 



[18] M. BOWNIK AND K. P. Ho, Atomic and molecular decomposition of anisotropic Triebel- 
Lizorkin spaces, Trans of Amer Math Soc, 385 (2005), pp. 1469-1510. 

[19] A. Calderon and A. Torchinsky, Parabolic maximal functions associted with a dis- 
tribution, Adv. Math., 24 (1977), pp. 101-171. 

[20] E. Candes and L. Demanet, The curvelet representation of wave propagators is opti- 
mally sparse. Comm. Pure Appl. Math., 58 (2005), pp. 1472- 1508. 

[21] M. Clausel and B. Vedel, Two optimality results about sample paths properties of 
Operator Scaling Gaussian Random Fields, Submitted, (2010). 

[22] S. Davies and p. Hall, Fractal analysis of surface roughness by using spatial data 
(with discussion)., J. Roy. Statist. Soc. Ser., B 61 (1999), pp. 3-37. 

[23] R. A. DeVore, S. V. Konyagin, and V. N. Temlyakov, Hyperbolic wavelet approx- 
imation, Constr. Approx., 14 (1998), pp. 1-26. 

[24] D. DONOHO, Wedgelets: Nearly minimax estimation of edges, Ann. Stat., 27 (1999), 
pp. 353-382. 

[25] K. Falconer, Fractal Geometry: Mathematical Foundations and Applications, Wiley 
and sons, 1990. 

[26] G. FoLLAND and E. Stein, Hardy spaces on homogeneous groups.. Mathematical Notes, 
28, Princeton University Press and University of Tokyo Press, 1982. 

[27] M. Frazier and B. Jawerth, Decomposition of besov spaces, Indiana University Math- 
ematics Journal, 34 (1985), pp. 777-799. 

[28] G. Garrigos and A. Tabacco, Wavelet decompositions of anisotropic Besov spaces. 
Math. Nachr., 239-240 (2002), pp. 80-102. 

[29] K. Guo and D. Labate, Representation of Fourier integral operators using shearlets. 
The Journal of Fourier Analysis and Applications, 14 (2008), pp. 327-370. 

[30] , Analysis and detection of surface discontinuities using the 3d continuous shearlet 

transform, Appl.Comp.Harm. Anal., 30 (2011), pp. 231-242. 

[31] D. Haroske and E. Tamasi, Wavelet frames for distributions in anisotropic Besov 
spaces, Georg. Math. J, 12 (2005), pp. 637-658. 

[32] R. Hochmuth, N-term approximation in anisotropic function spaces.. Math. Nachr., 
244 (2002), pp. 131-149. 

[33] , Wavelet characterizations for Anisotropic Besov Spaces, Appl. Comput. Harmon. 

Anal., 12 (2002), pp. 179-208. 

[34] L. Jacques, L. Duval, C. Chaux, and G. Peyre, A panorama on Multiscale Geo- 
metric Representations, Intertwining Spatial, Directional and Frequency Selectivity, Sub- 
mitted, (2012). 



34 



[35] S. Jaffard, Wavelet techniques in multifractal analysis, fractal geometry and applica- 
tions, Proc. Symp. Pure Math., 72 (2004), pp. 91-151. 

[36] , Pointwise and directional regularity of nonharmonic Fourier series, Appl. Comp. 

Harm. Anal., 28 (2010), pp. 251-266. 

[37] P. Lakhonchai, J. Sampo, and S. Sumetkijakan, Shearlet transforms and Holder 
regularities. International Journal of Wavelets, Multiresolution and Information Process- 
ing, To appear (2012). 

[38] T. Long and H. Triebel, Equivalent norms and Schauder bases in anisotropic besov 
spaces, Proc. Roy. Soc. Endinb., 84 A (1979), pp. 177-183. 

[39] Y. Meyer, Ondelettes et Operateurs : vol. I., Paris, Hermann, 1990. 

[40] M. Neumann, Multivariate wavelet thresolding in anisotropic function spaces, Statistica 
sinica, 10 (2000), pp. 399-431. 

[41] M. Neumann and R. von Sachs, Wavelet thresolding in anisotropic function classes 
and application to adaptative estimation of evolutionary spectra, Annals of Statistics, 25 
(1997), pp. 38-76. 

[42] K. Nualtong and S. Sumetkijakan, Analysis of Holder regularities by wavelet-like 
transforms with parabolic scaling,. Thai Journal of Mathematics, (2005), pp. 275-283. 

[43] G. Parisi and U. Frisch, On the singularity spectrum of fully developped turbulence, 
in Turbulence and Predictability in Geophysical Fluid Dynamics, Proceedings of the 
International Summer School in Physics Enrico Fermi,North Holland, (1985), pp. 84-87. 

[44] L. Ponson, D. Bonamy, H. Auradou, G. Mourot, S. Morel, E. Bouchaud, 
C. Guillot, and J. Hulin, Anisotropic self-affine properties of experimental fracture 
surfaces., Int.Journ. of fracture, 140 (2006), pp. 27-37. 

[45] S. Roux, M. Clausel, B. Vedel, S. Jaffard, and P. Abry, The Hyperbolic Wavelet 
Transform for self-similar anisotropic texture analysis. Submitted, (2012). 

[46] J. Sampo and S. Sumetkijakan, Estimations of Holder regularities and direction 
of singularity by hart smith and curvelet transforms. Journal of Fourier Analysis and 
Applications, 15 (2009), pp. 58-79. 

[47] H. Schmeisser and H. Triebel, Spaces of functions of mixed smoothness and their 
relations to approximation from hyperbolic crosses, Journal of Approximation Theory, 
128 (2004), pp. 115-150. 

[48] W. SiCKEL and T. Ullrich, Tensor products of Sobolev-Besov spaces and applications 
to approximation from the hyperbolic cross, Journal of Approximation Theory, 161 (2009), 
pp. 748-786. 

[49] H. Smith, A Hardy space for Fourier integral operators, J. Geom. Anal., 8 (1998), 
pp. 629-653. 

[50] C. Tricot, Two definitions of fractional dimension. Math. Proc. Cambridge Philos. 
Soc, 18 (1991), pp. 54-74. 



35 



[51] H. Triebel, Interpolation theory, function spaces, differential operators., Amsterdam, 
North-Holland, 1978. 

[52] , Wavelet bases in anisotropic function spaces., Function Spaces, Differential Opera- 
tors and Nonlinear Analysis, FSDONA-04. Milovy, Czech Republic, 18 (2004), pp. 529- 
550. 

[53] , Theory of function spaces III, Birhauser, Basel, 2006. 

[54] P. H. Westerink, Subband coding of images, PhD thesis. Delft University of Technol- 
ogy, Delft, The Netherlands, 1989. 



36 



